Detailed analysis of the gluonic excitation in the three-quark system in lattice QCD 



Tom T. Takahashi 

Yukawa Institute for Theoretical Physics, Kyoto University, 
Kitashirakawa-Oiwake, Sakyo, Kyoto 606-8502, Japan 



o 
o 

(N 

Or 
<D 



> 

in 
o 

o\ : 
o . 
^t- ■ 
o : 

<D . 
43 



X 



Hideo Suganuma 
Faculty of Science, Tokyo Institute of Technology, 
Ohokayama 2-12-1, Meguro, Tokyo 152-8551, Japan 
(Dated: February 1, 2008) 

We study the excited-state potential and the gluonic excitation in the static three-quark (3Q) 
system using SU(3) lattice QCD with 16 3 x 32 at /3=5.8 and 6.0 at the quenched level. For about 
100 different patterns of spatially-fixed 3Q systems, we accurately extract the excited-state potential 
Via' together with the ground-state potential V^ s ' by diagonalizing the QCD Hamiltonian in the 



presence of three quarks. The gluonic excitation energy Ai?3Q = Vac?' — V^' is found to be about 



3 q together with the ground-state potential V^q ' 

T '3Q ~ "3Q 

1 GeV at the typical hadronic scale. This large gluonic-excitation energy is conjectured to give a 
physical reason of the success of the quark model for low-lying hadrons even without explicit gluonic 
modes. We investigate the functional form of AEsq in terms of the 3Q location. The lattice data of 
AEsq are relatively well reproduced by the "inverse Mercedes Ansatz" with the "modified Y-type 
flux-tube length" , which indicates that the gluonic-excitation mode is realized as a complicated bulk 
excitation of the whole 3Q system. 



I. INTRODUCTION 

It is widely accepted that quantum chromodynamics 
(QCD) l] is the fundamental theory of the strong in- 
teraction for hadrons and nuclei. Nevertheless, it still 
remains as a difficult problem to derive low-energy phys- 
ical quantities directly from QCD in an analytic man- 
ner. The perturbative QCD calculation, which success- 
fully describes the high-energy process, cannot be applied 
at the hadronic scale, since the QCD coupling constant 
becomes large in the infrared region. Furthermore, the 
strong-coupling nature of QCD leads to a highly non- 
trival vacuum with rich nonperturbative phenomena such 
as color confinement and spontaneous chiral symmetry 
breaking. 

In this decade, the lattice-QCD Monte Carlo calcula- 
tion has been recognized as a reliable nonperturbative 
method for the quantitative analysis of QCD, and the 
nonperturbative analysis based on QCD has become one 
of the central and important issues in the hadron physics. 
For instance, lattice QCD calculations successfully repro- 
duce hadron mass spectra and also predict the QCD 
phase transition at finite temperatures and densities 
On the other hand, the underlying structure of hadrons 
is not yet well investigated using lattice QCD. In fact, 
the application of lattice QCD is just beginning for the 
research of the hadron structure and related excitation 
modes in terms of quarks and gluons. 

The investigation on the hadron structure and the ex- 
citation modes has a long history in the particle physics. 
In 1969, Y. Nambu first pointed out the string picture for 
hadrons 0, 0| to explain the Veneziano amplitude on 
hadron reactions and resonances. Since then, the string 
picture has been one of the most important scenarios for 
hadrons 0, and has provided many interesting ideas 
in the wide region of the elementary particle physics [t| . 



For instance, the hadronic string creates infinite num- 
ber of hadron resonances as the vibrational modes, and 
these excitations lead to the Hagedorn "ultimate" tem- 
perature 

[HHIl, which gives an interesting theoretical 
picture for the QCD phase transition. 

For real hadrons, the hadronic string has a spatial ex- 
tension like a "flux-tube" dEHH 0j , as the result of 
one-dimensional squeezing of the color-electric flux in ac- 
cordance with color confinement 0, 0] • Therefore, the 
vibrational modes of the hadronic flux-tube should be 
much more complicated, and the analysis of the excita- 
tion modes is important to clarify the underlying picture 
for real hadrons. 

In the language of QCD, such non-quark-origin excita- 
tion is called as the "gluonic excitation" [TH Ha Il7j , and 
is physically interpreted as the excitation of the gluon- 
field configuration in the presence of the quark-antiquark 
pair or the three quarks in a color-singlet state. 

In the hadron physics, the gluonic excitation is one of 
the interesting phenomena beyond the quark model, and 
relates to the hybrid hadrons such as qqG and qqqG in 
the valence picture. For instance, the hybrid meson in- 
cludes the exotic hadrons with the exotic quantum num- 
ber such as J PC = CP ~ , H ,1 '",2"' which cannot 
be constructed in the simple quark model ^^|. Then, it 
is important to investigate the gluonic excitation with 
lattice QCD, not only from the theoretical viewpoint but 
also from the experimental viewpoint. 

In this paper, we study the excited-state three-quark 
(3Q) potential and the gluonic excitation in baryons us- 
ing lattice QCD [15|, to get deeper insight on these exci- 
tations beyond the hypothetical models such as the string 
and the flux-tube models. In QCD, the excited-state 3Q 
potential is defined as the energy of the excited state 
of the gluon-field configuration in the presence of the 
static three quarks, and the gluonic-excitation energy is 
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expressed as the e nergy difference between the ground- 
state 3Q potential [lO 120 |21| and the excited-state 3Q 
potential. 

Note that the inter-quark potential is one of the most 
fundamental quantities directly connected to color con- 
finement, and plays the important role for hadron prop- 
erties. As for the Q-Q system, a lot of lattice QCD stud- 
ies [2E HE 13 HH have shown that the ground-state Q-Q 
potential Vqq is well described by the Coulomb plus lin- 
ear potential, VJ?J = — — t^+Cqq^+Cqq, with the inter- 
quark distance r. Its behavior at short distances can 
be explained by the Coulomb interaction from the one- 
gluon-exchange (OGE) process, and the linear confine- 
ment term seems to indicate the flux-tube picture 0,0], 
where the quark and the antiquark are linked by one di- 
mensional flux-tube with the string tension of a ~ 0.89 
GeV/fm. 

In QCD, the three-body force among three quarks is 
also a "primary" force reflecting the SU(3) gauge symme- 
try, while the three-body force is regarded as a residual 
interaction in most fields in physics. In fact, the 3Q po- 
tential is directly responsible for the structure and prop- 
erties of baryons, similar to the relevant role of the Q-Q 
potential for meson properties, and both the Q-Q poten- 
tial and the 3Q potential are equally important funda- 
mental quantities in QCD. Furthermore, the 3Q poten- 
tial is the key quantity to clarify the quark confinement 
mechanism in baryons. However, in contrast with the 
Q-Q potential, there were only a few preliminary lattice- 
QCD works [M El HI done in 1980's for the ground- 
state 3Q potential before our first study in 1999 |19| . 

Since 1999, we have performed accurate calculations 
and detailed analyses of the ground-state 3Q poten- 
tial V^q' in lattice QCD with the smearing method for 
more than 300 different patterns of 3Q systems |2l|. In 

Refs. minim, 

we have shown that the ground-state 
3Q potential V^q' is well described by the Coulomb plus 
Y-type linear potential, i.e., the Y-Ansatz, 
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within 1%-level deviation. Here, L m i n denotes the min- 
imal value of the total leng th of color-flux-tubes linking 
the three quarks [ll 111 120121 129ll30j. which is schemat- 
ically illustrated in Fig. We have found two remark- 
able feature's, the universality of the string tension as 
ct 3 q ~ Oqq and the OGE result as A 3 q ~ ^qq- (Very 
recently, we show that the multi-quark potential is well 
described by the OGE Coulomb plus multi-Y type linear 
potential |3l|, which also supports the Y-Ansatz.) 

We briefly summarize several other recent studies on 
the ground-state 3Q potential to clarify the current sta- 
tus of the Y-Ansatz. 

After our first study, de Forcrand's group tested only 
about 20 equilateral 3Q configurations in lattice QCD, 
and supported the A-Ansatz [33]. However, their data 
seem to indicate the Y-Ansatz with an overall constant 



Lmin=PQl+PQ2+PQ3 



FIG. 1: The Y-type flux-tube configuration of the 3Q sys- 
tem with the minimal value of the total flux-tube length, 
Lmin — X^i-i PQi- There appears a physical junction linking 
the three flux tubes at the Fermat point P. 



shift, and the deviation seems to originate from the data 
at very short distances, where the linear potential is neg- 
ligible compared with the Coulomb contribution. (Also, 
their usage of the continuum Coulomb potential may be 
problematic for the lattice data at very short distances.) 
Recently, de Forcrand's group seems to c hang e their opin- 
ion from the A-Ansatz to the Y-Ansatz [33| . 

One of the theoretical basis of the A-Ansatz was Corn- 
wall's conjecture based on the vortex vacuum model [34| . 
Very recently, motivated by our studies, Cornwall re- 
examined his previous work and found an error in the 
model calculation. The corrected answer was the Y- 
Ansatz instead of the A-Ansatz 35] . 

As another analytical work, Kuzmenko and Simonov 
also showed that the Delta-shape is impossible from 
gauge-invariance point of view, and the Y-shaped con- 
figuration is the only possible for the 3Q system |3^ . 

As a clear evidence for the Y-Ansatz, Ichie et al. per- 
formed the direct measurement on the action density in 
the spatially-fixed 3Q system in lattice QCD using the 
maximally- Abelian projection, and observed a clear Y- 
type flux-tube profile |16ll37| . 

In this way, the Y-Ansatz is also supported by various 
studies of other groups [331 HE HE 113 ■ and therefore 
the Y-Ansatz for the ground-state 3Q potential is almost 
settled both in lattice QCD and in analytic framework. 

As for the excited-state 3Q potential, however, there is 
no lattice QCD study besides our previous work [l5|. In 
this paper, we present the detailed analysis of the excited- 
state 3Q potential and the gluonic excitation for about 
100 different patterns of the 3Q static systems in SU(3) 
lattice QCD at the quenched level. 

This paper is organized as follows. In section [HI we 
give a necessary formalism on the relation between the 
3Q Wilson loop and the QCD Hamiltonian. We then give 
the lattice QCD formalism to extract the excited-state 
3Q potential in section IIIII and show the lattice QCD 
results in section IIVI In section we investigate the 
functional form of the gluonic excitation energy from the 
lattice QCD data. In section lv*Tl we discuss the physical 
implication of the obtained lattice results, and consider 
the physical reason of the success of the quark model for 
low-lying hadrons in terms of the gluonic excitation. Sec- 
tion is devoted to the summary and the conclusion. 
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II. 3Q WILSON LOOP AND 3Q POTENTIALS 
A. 3Q Wilson loop and QCD Hamiltonian 

To begin with, let us consider the physical modes in 
the static 3Q system. Wc denote the nth excited state 
by \n) (n=0,l,2,3,..) for the physical eigenstates of the 
QCD Hamiltonian H for the spatially-fixed 3Q system. 
For the simple notation, the ground state is expressed 
as the "Oth excited state" in this paper. Since the three 
quarks are spatially fixed in this case, the eigenvalue of 
H is expressed by the static 3Q potential as 



H\n) =E n \n) = V n \n), 



(2) 



where V n denotes the nth excited-state 3Q potential. We 
take the normalization condition as (m\n) — S mn . Note 
that both V n and \n) are universal physical quantities 
related to the QCD Hamiltonian H. 




FIG. 2: The 3Q Wilson loop W 3Q . The gauge-covariant 3Q 
state |$) is generated at t = and is annihilated at t = T. 
The three quarks are spatially fixed in R 3 for < t < T. 



Similar to the calculation of the Q-Q potential with 
the Wilson loop, the 3Q potential can be calculated with 
the 3Q Wilson loop W3Q defined as 



W; 



3Q 
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(3) 



with U k = Pexp{ig / r dx^A^x)} (k = 1, 2, 3). Here, P 
denotes the path-ordered product along the path denoted 
by T fc in Fig.lU 

In the 3Q Wilson loop, a gauge-covariant 3Q state |$) 
is generated at t = and is annihilated at t = T with the 
three quarks being spatially fixed in R 3 for < t < T. In 
general, the 3Q state |$) in the 3Q Wilson loop is not an 
eigenstate of the QCD Hamiltonian H for the spatially- 
fixed 3Q system, and can be expressed with a linear com- 
bination of the 3Q physical eigenstates \n) (n=0,l,2..) as 



|$) =c |0) + Cl |l)+c 2 |2) + 



(4) 



where the complex coefficients c n satisfy X^^Li l c «| 2 = 1 
and express the overlap with each state \n) in the 3Q 
Wilson loop. 



Since the Euclidean time evolution of the 3Q state 
|$(£)} is expressed with the operator e~ Ht , which cor- 
responds to the transfer matrix in lattice QCD, the ex- 
pectation value of the 3Q Wilson loop is expressed as 



(W 3Q (T)) = ($(T)|$(0)) = <$| e -" J |$) 



(5) 
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) = E 



c n ?e~ v - T . 



n=0 



Then, the low- lying potentials, e.g., the ground-state 
potential V^' = Vq and the 1st excited-state potential 
VJjq ' = V\ , can be obtained from the large-T behavior of 
(W3q(T)) where higher excited-state contributions are 
negligible. In practical lattice simulations, however, it 
is difficult to calculate (W3q(T)) accurately for large T, 
since its value decreases exponentially with T. 

Therefore, for the accurate calculation of low-lying po- 
tentials, it is desired to reduce higher excited-state com- 
ponents in the 3Q state |$) in the 3Q Wilson loop. 



B. The smearing method 

For this purpose, we use the smearing method |3S| to 
reduce highly excited-state components in the 3Q state 
$). The smearing is just a method for the state com- 
position in a gauge-covariant manner, and never changes 
physical quantities and gauge configurations, unlike the 
cooling. 

The smearing for link-variables is expressed as the 
iterative replacement of the spatial link- variable Ui(s) 
(i = 1,2,3) by the obscured link-variable Ui(s) e SU(3) 
which maximizes 



aUi(s) + Y,{U 3 (s)U t {s+j)u}(s + i) 



Re trj Uj(s) 

+U}(s-j)U i (s-j)U j (s + i-j)} 



(6) 



leaving the temporal link- variable K^s) unchanged. 

Note that the smeared 3Q Wilson loop composed by 
the smeared link-variable is expressed as a spatially- 
extended operator in terms of the original link-variable 
[/^(s), and therefore the 3Q state |$) and the coefficients 
c n in the smeared 3Q Wilson loop are changed according 
to the iteration number N sml of the smearing. In other 
words, the coefficients c„ in the smeared 3Q Wilson loop 
can be controlled to some extent in the smearing pro- 
cedure, by changing the iteration number -/V smr and the 
smearing parameter a. 

For instance, as the iteration number -/V smr of the 
smearing increases from iV smr =0, the excited-state com- 
ponents in the iV smr th smeared 3Q Wilson loop gradually 
decrease, and finally the 22th (42th) smeared 3Q Wilson 
loop is almost ground-state saturated as[$) ~ c§ 1 0} for 
a = 2.3 on the lattice with 0=5.8 (6.0) |2j. Then, with 
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the properly smeared 3Q Wilson loop (Wsq(T)), V^q' 
can be accurately calculated as V^q' ~ — ^ ln(W / 3 Q(T)) 
even at relatively small T [20L [2l| . 

The smearing method is also useful to extract the low- 
lying 3Q potentials, i.e., the ground-state potential V^q S ' 
and the 1st excited-state potential Vgq' [l5llTfil |. 

III. FORMALISM 

In this section, we present the formalism of the vari- 
ational method [15J, |3JJ |4CJ to extract the excited-state 
potential and the gluonic excitation by diagonalizing the 
QCD Hamiltonian in the presence of static quarks. 

Suppose that |$fe) (k = 0, 1, 2, 3, ..) are arbitrary given 
independent 3Q states for the spatially-fixed 3Q sys- 
tem. In general, each 3Q state can be expressed 
with a linear combination of the 3Q physical eigenstates 
\n) (n = 0, 1,2, ..) as 

\<P k )=c k \0) + c k 1 \l)+c k 2 \2) + .. (7) 

Here, the coefficients c k depend on the selection of the 3Q 
states \&k), and hence they are not universal quantities. 
(Unlike the Q-Q system, there is no definite symmetry in 
the 3Q system, so that we do not construct the 3Q states 
which carry the specific quantum number.) 

The Euclidean time-evolution of the 3Q state 
is expressed with the operator e~ Ht , which corresponds 
to the transfer matrix in lattice QCD. The overlap 
(3?j(T)|<I>fc(0)) is given by the generalized 3Q Wilson loop 
sandwiched by the initial state at t = and the final 
state <frj at t = T, and is expressed as the correlation 
matrix in the Euclidean Heisenberg picture: 

wf = ($ J -(T)|$ fc (o)) = ($ j |w 3 Q(r)|$ ifc ) 

00 00 

= (^\e-^ k )=J2J24 n c k n (rn\e- 6T \n) 

m—0 n— 

oc 

= J2? n c k e- V ~ T . (8) 

We define the matrix C and the diagonal matrix At by 

C nk = c k , A T m = e- v " T 6 mn , (9) 

and rewrite the above relation as 

W T = CUtC. (10) 

Note here that C is not a unitary matrix, and therefore 
this relation does not mean the simple diagonalization 
by the unitary transformation, which connects the two 
matrices by the similarity transformation. 

Since we are interested in the 3Q potential V n (n = 
0, 1, 2, ..) in At rather than the non-universal matrix C, 
we single out V n from the 3Q Wilson loop Wt using the 



following prescription. From Ea. (|10fl . we obtain 

W^Wt+i = {C^AtC^C^At+iC 
= C^A^At+iC 
= C- 1 disig(e- Vo ,e- v \..)C. (11) 

Now, Wt+x is expressed by the similarity transfor- 
mation of the diagonal matrix diag(e~ Vo , e~ Vl , ..). There- 
fore, e~ Vn can be obtained as the eigenvalues of the ma- 
trix Wt+i, i.e., the solutions of the secular equation, 

detiW^Wr+i - tl} = det{A 7 ; 1 A T+1 - tl} 

= H(e- v --t) = 0, (12) 

n 

where / denotes the unit matrix. The largest eigenvalue 
corresponds to e~ Vo and the nth largest eigenvalue cor- 
responds to e~ Vn . 

In this way, the 3Q potential V n (n — 0,1,2,..) can 
be obtained from the matrix W t ~ 1 Wt+i- In the practi- 
cal calculation, we prepare N independent sample states 
|$fc) (k — 0, 1, ..,N — 1). If one chooses appropriate 
states |$fc) which does not include highly excited-state 
components, one can truncate the physical states as 
\n) (n = 0, 1,2, ..,JV - 1). Then, W T , C and A T are 
reduced into N x N matrices, and the secular equation 
Ea. (|12|) becomes the iVth order equation. 

IV. LATTICE QCD RESULT 

In this section, we show the lattice QCD result of the 
1st excited-state 3Q potential U 3 e q ■ and the gluonic exci- 
tation energy AE^q as well as the ground-state potential 
U 3 q S ' for the spatially-fixed static 3Q system. The SU(3) 
lattice QCD calculation is done with the standard plaque- 
tte action on 16 3 x 32 at (3—5.8 and 6.0 at the quenched 
level. The lattice spacing is found to be a ~ 0.15 fm at 
j3 = 5.8 and a ~ 0.1 fm at j3 = 6.0, which are set to 
reproduce the string tension (7=0.89 GeV/fm in the Q-Q 
potential In Table [Q we summarize the simulation 
condition and related information of the present lattice 
QCD calculation for the 3Q potentials. 

TABLE I: The simulation condition and related information. 
For each 0, we list the corresponding lattice spacing a, the 
lattice size, the number A3Q of the different patterns of the 3Q 
system analyzed, the gauge configuration number N con f used 
for the measurement, the number Atherm of sweeps for the 
thermalization, the number A aep of sweeps for the separation, 
the smearing parameter a and the iteration number A amr used 
for the extraction of the 3Q potentials, U|q S ' and V£q' ■ 

P a [fin] lattice 

5.8 0.15 16 3 x 32 24 200 10,000 500 2.3 8,12,16,20 
6.0 0.10 16 3 x 32 73 149 10,000 500 2.3 16,24,32,40 



From now on, we concentrate ourselves on the ground 
state |0) and the 1st excited state |1) in the spatially-fixed 
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3Q system. To extract Vq and V±, we need to prepare at 
least two independent states \&k)(k = 0,1), and con- 
struct the 2x2 matrix W^ 1 Wt+i with them. Here, the 
sample states \&k) can be freely chosen, as long as they 
satisfy the two conditions: the linear independence and 
the smallness of the higher excited-state components \n) 
with n > 2, which leads to |$ fe ) ~ cg|0) + cf |1). 

As the sample 3Q states $fc), we adopt the properly 
smeared 3Q states since the higher excited-state com- 
ponents are reduced in them [2l|. Here, the smearing 
parameter is fixed to be a — 2.3. After some numeri- 
cal check on the above two conditions, we adopt the 8th, 
12th, 16th, 20th smeared 3Q states at /3 = 5.8 and the 
16th, 24th, 32nd, 40th smeared 3Q states at (3 = 6.0 as 
the candidates of the sample 3Q states. Owing to the 
intervals of 4 (8) iterations at (3=5.8 (6.0), these smeared 
states are clearly independent of each other. The iV smr th 
smeared state with N smi > 8 (16) at (3=5.8 (6.0) has 
small higher excited-state components. 
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FIG. 3: An example of the effective-mass plot for Vo(T) (open 
circles) and Vi(T) (filled circles) obtained from the eigen- 
values of W^ 1 Wt+i at each T for the spatially-fixed three 
quarks put on (1, 0, 0), (0, 1, 0) and (0, 0, 1) in the lattice unit 
at (3=5.8. For V {T) and Vi(T) at each T, we plot 6 data 
obtained from the 6 pairs of the state combination. 

For each possible pairing (j, k) of these 4 states, 
we calculate the generalized 3Q Wilson loop Wj- = 
($ J (T)|$ fc (0)), and evaluate V and Vx with Eq.JT^I. We 
plot in Fig. |3 an example of the "effective mass" plot for 
V (T) and Vi(T) obtained from Eq.iO at (3 = 5.8 as the 
function of the temporal separation T. For the estima- 
tion of the statistical error of the lattice data, we adopt 
the jack-knife error estimate. 

For V (T) and Vi(T) at each T, we have 6 (= 4 C 2 ) data 
obtained from 6 pairs of the combination among the 4 
sample states, i.e., the 8th, 12th, 16th, 20th smeared 3Q 
states. As shown in Fig. |3J these 6 data almost coincide 
and the T-dependence of Vo (T) and V\ (T) is rather small 
in a certain region of T. This indicates the smallness of 
the higher excited-state components |n) with n > 2 in the 
sample states, since such contaminations lead a nontrivial 
T-dependence in Vq (T) and V\ (T) and make the stability 
lost. 

For the accurate measurement, we select the best 



pairing providing the most stable effective-mass plot, 
which physically means the smallest contamination of the 
higher excited states in them. (If the effective-mass plot 
does not show a plateau, we exclude the 3Q configura- 
tion from the analysis to keep the accuracy.) With the 
selected two states, we extract the ground-state poten- 
tial V^q' and the 1st excited-state potential V^q ' using 
the x 2 At as Vf Q s - = V Q {T) and V££ = Vi(T) in the fit 
range of T, where the plateau is observed. We perform 
the above procedure for each 3Q system. 

In TableslHllHllandllVl we summarize the lattice QCD 
results for the ground-state 3Q potential V^q' and the 1st 
excited-state potential V^ e q ■ at the quenched level. Here, 
we analyze the 97 different patterns of the spatially-fixed 
3Q system in total on the 16 3 x 32 lattice at (3 =5.8 and 
6.0. 

TABLE II: The ground-state 3Q potential VJq 8, and the 1st 
excited-state 3Q potential VjJq ' at (3 = 5.8 in the lattice unit. 
The label (l,m,n) denotes the 3Q system where the three 
quarks are put on (la, 0,0), (0,ma,0) and (0,0, na) in R 3 . 



I 1 TY1 T) \ 

I I ^ / / b , lb) 


*3Q 


T/g.S. 

V 3Q 


AE :iQ = V£& - V 3 % s - 


(0,1,1) 


1.9816( 95) 


0.7711( 3) 


1.2104 


95) 


(0,1,2) 


1.9943( 72) 


0.9682( 4) 


1.0261 


72) 


(0,1,3) 


2.0252( 92) 


1.1134( 7) 


0.9118 


90) 


(0,2,2) 


2.0980( 80) 


1.1377( 6) 


0.9603 


81) 


(0,2,3) 


2.1551( 87) 


1.2686( 9) 


0.8866 


86) 


(0,3,3) 


2.2125(114) 


1.3914(13) 


0.8211( 


112) 


(1,1,1) 


2.0488( 90) 


0.9176( 4) 


1.1312 


90) 


(1,1,2) 


2.0727( 75) 


1.0686( 5) 


1.0041 


75) 


(1,1,3) 


2.1023( 73) 


1.2004( 7) 


0.9019 


73) 


(1,1,4) 


2.1580( 93) 


1.3201(10) 


0.8380 


92) 


(1,2,2) 


2.1405( 72) 


1.1907( 7) 


0.9498 


71) 


(1,2,3) 


2.1899( 71) 


1.3084( 9) 


0.8815 


70) 


(1,2,4) 


2.2516( 79) 


1.4221(12) 


0.8296 


78) 


(1,3,4) 


2.2907( 91) 


1.5260(15) 


0.7647 


88) 


(1,4,4) 


2.3807(138) 


1.6322(20) 


0.7485( 


136) 


(2,2,2) 


2.1776(111) 


1.2844(10) 


0.8932( 


110) 


(2,2,3) 


2.2242( 96) 


1.3882(11) 


0.8360 


95) 


(2,2,4) 


2.2799( 98) 


1.4952(15) 


0.7847 


99) 


(2,3,4) 


2.3637(100) 


1.5853(18) 


0.7784 


99) 


(2,4,4) 


2.4108(137) 


1.6836(23) 


0.7271( 


135) 


(3,3,3) 


2.3408(168) 


1.5680(19) 


0.7728( 


166) 


(3,3,4) 


2.3958(151) 


1.6635(22) 


0.7323( 


146) 


(3,4,4) 


2.4645(177) 


1.7565(30) 


0.7081( 


173) 


(4,4,4) 


2.5245(340) 


1.8408(42) 


0.6837( 


343) 



The lattice data of Vf Q s - in Tables |H| |IH] and |IV| also 
indicate the validity of the Y-Ansatz for the ground-state 
3Q potential. Note that the present data of V^ 8 ' are con- 
sidered to include almost no excited-state contributions, 
since they are extracted by diagonalizing the correlation 
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TABLE III: The ground-state 3Q potential V 3 g Q s - and the 1st 
excited-state 3Q potential V 3 c q ' at /3 = 6.0 in the lattice unit. 
The label (Z, m, n) denotes the 3Q system where the three 
quarks are put on (la, 0,0), (0, ma,0) and (0,0, na) in R' 5 . 



!, 771, 71) 




^3Q 












_ T/G.S. 


0,1,1 ) 


1 


5973 ( 


701) 





6765 


; 6) 





9209( 


702) 


0,1,2 ) 


1 


6502( 


190) 





8233 


; 6) 





8269( 


190) 


0,1,3 ) 


1 


6566 


78) 





9159 


: 6) 





7407 


77) 


0,1,4 ) 


1 


6762( 


101) 





9868 


; 9) 





6893( 


100) 


0,1,5 ) 


1 


6861 


92) 


1 


0491 


12) 





6370 


91) 


0,1,6 ) 


1 


7135 


99) 


1 


1062 


17) 





6073 


98) 


0,2,2 ) 


1 


7380 


75) 





9452 


: 6) 





7929 


74) 


0,2,3 ) 


1 


7559 


79) 


1 


0266 


: 7) 





7292 


78) 


0,2,5 ) 


1 


7858 


88) 


1 


1525 


13) 





6333 


87) 


0,2,6 ) 


1 


8140 


91) 


1 


2084 


16) 





6055 


90) 


0,3,3 ) 


1 


7880 


65) 


1 


1001 


; 9) 





6879 


64) 


0,3,4 ) 


1 


8108 


70) 


1 


1620 


ii) 





6488 


70) 


0,3,5 ) 


1 


8325 


78) 


1 


2191 


15) 





6134 


77) 


0,3,6 ) 


1 


8631 


83) 


1 


2721 


18) 





5910 


83) 


0,4,4 ) 


1 


8521 


75) 


1 


2198 


14) 





6323 


74) 


0,4,5 ) 


1 


8843 


77) 


1 


2742 


17) 





6101 


76) 


0,4,6 ) 


1 


9074 


92) 


1 


3281 


20) 





5793 


91) 


0,5,5 ) 


1 


9093 


76) 


1 


3283 


19) 





5809 


76) 


0,5,6 ) 


1 


9450 


84) 


1 


3791 


24) 





5659 


84) 


0,6,6 ) 


1 


9766 ( 


100) 


1 


4286 


25) 





5480 


98) 


1,1,1 ) 


1 


7035 


66) 





7941 


[ 3) 





9094 


65) 


1,1,2 ) 


1 


7236 


77) 





8994 


' 4) 





8242 


76) 


1,1,3 ) 


1 


7196 


74) 





9817 


I 6) 





7378 


73) 


1,1,4 ) 


1 


7348 


97) 


1 


0494 


s 9) 





6854 


96) 


1,1,5 ) 


1 


7422 


84) 


1 


1105 


12) 





6317 


83) 


1,2,2 ) 


1 


7476 


58) 





9810 


' 5) 





7666 


58) 


1,2,3 ) 


1 


7693 


67) 


1 


0521 


v 7 ) 





7171 


66) 


1,2,4 ) 


1 


7856 


73) 


1 


1151 


10) 





6705 


71) 


1,2,5 ) 


1 


8031 


80) 


1 


1741 


13) 





6290 


78) 


1,2,6 ) 


1 


8273 


87) 


1 


2279 


16) 





5994 


85) 


1,3,3 ) 


1 


7964 


59) 


1 


1161 


19) 





6804 


58) 


1,3,4 ) 


1 


8213 


66) 


1 


1745 


11) 





6467 


64) 


1,3,5 ) 


1 


8123( 


202) 


1 


2286 


22) 





5837( 


202) 


1,3,6 ) 


1 


8800 


83) 


1 


2842 


19) 





5958 


82) 


1,4,4 ) 


1 


8483 


66) 


1 


2288 


14) 





6196 


64) 


1,4,5 ) 


1 


8882 


66) 


1 


2829 


18) 





6053 


64) 



TABLE IV: The ground-state 3Q potential and the 1st 
excited-state 3Q potential V^q' at /3 = 6.0 in the lattice unit. 
The notations are the same as in Table UTTl 



(Z, m, n) 


V£q 




T/g.S. 


AE 3Q = Vi£ - V£ Q B - 


(1,5,5 ) 


1.9210 


72) 

1 u J 


1.3343(19) 


d c Oct 


71) 


(1,5,6 ) 


1.9460 


{ 85) 


1.3843(23) 


0.5617 


' 84) 


(1,6,6 ) 


1.9855 


' 88 ) 


1.4328(26) 


0.5527 


[ 84) 


(2,2,2 ) 


1.7687 


' 61) 


1.0392( 6) 


U. (29t> 


DO J 


(2,2,3 ) 


1.7901 


; 58) 


1.0993( 8) 


0.6908 


, 57) 


(2,2,4 ) 


1.8107 


' 64) 


1.1575(10) 


0.6532 


06) 


(2,2,5 ) 


1.8290 


' 73) 


1.2129(13) 


0.6161 


71) 


(2,2,6 ) 


1.8631 


' 76) 


1.2667(17) 


0.5964 


74) 


(2,3,3 ) 


1.8109 


' 56 ) 


1.1509( 9) 


U.6600 


' C A \ 

04J 


(2,3,4 ) 


1.8354 


' 65) 


1.2044(13) 


0.6310 


RA\ 

64) 


(2,3,5 ) 


1.8632 


' 74) 


1.2575(16) 


0.6058 


; 73) 


(2,3,6 ) 


1.9053 


' 77) 


1.3096(18) 


U.o9o I 


74) 


(2,4,4 ) 


1.8575 


' 63) 


1.2547(15) 


U.6028 


62 ) 


(2,4,5 ) 


1.9045 


' 66) 


1.3068(18) 


U.o9r 1 


' a a \ 
64 J 


(2,4,6 ) 


1.9167( 


298) 


1.3542(34) 


0.5625 


297) 


(2,5,5 ) 


1.9295( 


300) 


1.3509(31) 


0.5786( 


301) 


(2,5,6 ) 


1.9689 


' 76) 


1.4037(24) 


U.o6o3 


' ta ^ 
74) 


(2,6,6 ) 


1.9813 


' 85) 


1.4493(27) 


0.5320 


80) 


(3,3,3 ) 


1.8434 


' 55 ) 


1.1968(12) 


U.64D6 


06) 


(3,3,4 ) 


1.8695 


' 60) 


1.2467(14) 


U.6228 


59) 


(3,3,5 ) 


1.8923 


' 66 ) 


1.2963(18) 


0.5961 


[ 63) 


(3,3,6 ) 


1.9371 


; 69) 


1.3479(20) 


U.089z 


01 ) 


(3,4,4 ) 


1.8464( 


244) 


1.2879(24) 


U.0084 


O A A \ 

244 ) 


(3,4,5 ) 


1.9164 


' 67) 


1.3380(20) 


0.5784 


65) 


(3,4,6 ) 


1.9389 


' 77) 


1.3881(24) 


0.5507 


' 71 \ 

71) 


(3,5,5 ) 


1.9032( 


322) 


1.3794(34) 


U.0238 


o i n\ 
319J 


(3,5,6 ) 


1.9691 


78) 


1.4314(25) 


0.5377 


; 74) 


(3,6,6 ) 


1.9743( 


383) 


1.4717(51) 


0.5026( 


381) 


(4 4.4 ) 


1.9215 


' 61) 


1.3322(18) 


0.5893 


; 58) 


(4,4,5 ) 


1.9321 


73) 


1.3766(23) 


0.5555 


; 72) 


(4,4,6 ) 


1.9848 


' 72) 


1.4253(24) 


0.5595 


;7i) 


(4,5,5 ) 


1.9569( 


383) 


1.4190(44) 


0.5379( 


382) 


(4,5,6 ) 


1.9749( 


398) 


1.4582(54) 


0.5166( 


401) 


(4,6,6 ) 


1.9667( 


441) 


1.5039(65) 


0.4628( 


431) 


(5,5,6 ) 


2.0076( 


542) 


1.5050(61) 


0.5026( 


542) 


(5,6,6 ) 


2.0150( 


563) 


1.5445(63) 


0.4705( 


556) 


(6,6,6 ) 


2.0970( 


795) 


1.5925(79) 


0.5046( 


780) 



matrix ($j|e — |<&fc) in terms of the physical basis \n). 
As a consequence, the accuracy of the present data is 
better than that of the data in Refs. pol l2l|. 

In Fig. 0] we plot the ground-state potential V^q' and 
the 1st excited-state potential V 3 e q' as the function of 
the minimal length L m m of the Y-type flux-tube in the 
lattice unit. The open symbols are for the ground-state 
potential V^f', and the filled symbols are for the excited- 



state potential VJq in the 3Q system. (Note here that we 
use Lmin simply for the distinction between the different 
3Q configurations, and then symbols are not necessarily 
required to lie on a single curve in the figures.) Here, V 3 ^ s ' 
and V^q " are the lowest and the next-lowest eigenvalues 

of the QCD Hamiltonian H for the static 3Q system, 
and correspond to the ground-state and the 1st excited- 
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state energies induced by three static quarks in a color- 
singlet state. We note that the lattice results at /3 = 5.8 
and P = 6.0 well coincide in the physical unit besides an 
irrelevant overall constant. The gluonic excitation energy 
is expressed as AE 3 q 
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FIG. 5: The lattice QCD results of the gluonic excitation 



energy AE^q = V^q' — V£q ' plotted against L m i n , the minimal 
total length of flux-tubes linking the three quarks. The results 
at /3 = 5.8 and /3 = 6.0 well coincide in the physical unit. 



FIG. 4: The lattice QCD results of the ground-state 3Q po- 
tential V^q' (open circles) and the 1st excited-state 3Q po- 
tential V 3 c q ' (filled circles) plotted against Z/ m i n , the minimal 
total length of flux-tubes linking the three quarks. The lattice 
results at /3 = 5.8 and (3 = 6.0 well coincide in the physical 
unit besides an irrelevant overall constant. 

It is worth mentioning that the absolute values of the 
potentials cannot be determined in lattice QCD without 
ambiguity. In fact, the energy of the 3Q system measured 
with the 3Q Wilson loop contains an irrelevant constant 
term C3Q, which corresponds to the self-energies of the 
three static quarks under the lattice cutoff a -1 and di- 
verges in the continuum limit as a — > |2£|[2l|. However, 
the energy gap between any pair of two states does not 
suffer from the ambiguity and has physical meaning. In 
particular, the energy gap AE 3 q = V£q — V 3 S q' between 
the ground-state and excited-state potentials has definite 
physical meaning as the lowest gluonic excitation energy, 
and can be determined in lattice QCD without the am- 
biguity. 

Figure [5] shows the gluonic excitation energy AE 3 q = 
V£q ■ — V^q' as the function of L m i n . As a nontrivial fact, 
AE3Q is almost reproduced with a single-valued func- 
tion of Z/ m in, the minimal total length of the flux-tube 



in the 3Q system. This implies that the gluonic excita- 
tion energy AE$q is controlled by the whole size of the 
3Q system. This will be discussed in detail in the next 
section. 



V. FUNCTIONAL FORM OF GLUONIC 
EXCITATION ENERGY 

In this section, we investigate the functional form of the 
gluonic excitation energy Al?3Q(ri, t^, V3) in the static 
3Q system in terms of the 3Q location (z=l,2,3) using 
the lattice QCD data for AE 3 q at ,9=5.8 and 6.0. 

The 3Q static potential V3Q generally depends on the 
three independent variables which indicate the 3Q trian- 
gle, e.g., {a,b,c} for the three sides of the 3Q triangle, 
while the Q-Q potential V^q depends only on the rela- 
tive distance r. Therefore, the search for the functional 
form of V 3 c q' or AE$q is much more difficult than in the 
Q-Q case. 

Furthermore, unlike the ground-state 3Q potential 
VgQ S ', there are no clear theoretical candidates for the 
functional form of the excited-state 3Q potential V30 ' or 
the gluonic excitation energy AE^q. Hence, it is rather 



difficult to specify their functional forms, and we are 
obliged to perform "trial and error". Here, we consider 
various possible functional forms, and perform the x 2_n ^ 
for the lattice QCD data for each form. 

Since the Coulomb part originated from the OGE 
process is expected to be equal in V£q and V 3 q', the 
Coulomb contribution is considered to be cancelled in 
the combination of AE^q ee V 3 c q' — VJq 5 '. Then, the 
functional form of AE^q is expected to be simpler than 
the excited-state 3Q potential V£q'. Therefore, we in- 
vestigate the gluonic excitation energy AE 3 q in detail 
instead of VJj^' . 



A. Comparison with the Q-Q system 

To begin with, we attempt to express the gluonic 
excitation energy AE 3 q of the 3Q system in terms of 
Ai?qq(r) of the Q-Q system, with considering the physi- 
cal structure of the 3Q system. Here, we refer to Ref. 
on the lattice QCD data of A_Eq>q(r) of the Q-Q system. 

Let us consider the 3Q system as shown in Fig. ^ with 
the quark location Qi (z=l,2,3) and the Fermat point P. 
We denote the lengths of the three sides by a, b and c. 

If the Y-junction exhibits the fixed-edge nature, the 
gluonic excitation of the 3Q system is expected to re- 
semble that of the Q-Q system, since the three static 
quarks also play the role of the fixed edges. If it is the 
case, the lowest gluonic excitation energy AE 3 q in the 
3Q system would be expressed as 

AE 3Q ~ AE QQ (r = maxCPQ^PQ^PQa)). (13) 

If the excitation mode can be expressed as a vibrational 
mode on one side of the 3Q triangle, the lowest gluonic 
excitation is expressed by the lowest vibrational mode of 
the Q-Q flux as 

AE 3Q ~ AE QQ (r = max(a, b, c)). (14) 

However, these fits cannot reproduce the lattice QCD 
data of A-E3Q at all. Then, the 3Q gluonic excitation is 
considered as the bulk excitation of the whole 3Q system, 
rather than the excitation of its partial system. 

We have also checked many trial forms with AEqq (r) 
such as 

AE 3Q ~ AE QQ (r = L min ), (15) 
AE 3Q ~ AE QQ (r = a + b + c). (16) 

However, all of them fail to reproduce the lattice QCD 
result of AE 3 q. 



B. The inverse Mercedes Ansatz 

Next, as a trial, we attempt to plot the gluonic ex- 
citation AE3Q in the 3Q system against the minimal 



S 



Q1 




FIG. 6: The "Mercedes form" for the 3Q system, where 
quarks are put on Qi (i=l,2,3). The modified Y-type flux- 
tube length L Y is expressed as L Y = 4 T^j-j-j P»Q.7- P de- 
notes the Fermat point of AQ1Q2Q3, and Pi is taken on the 
extended line of QiP so as to satisfy PPi = £. The "inverse 
Mercedes Ansatz" is defined as A.E3Q = K/L Y + G. 

total length of the flux-tubes linking the three quarks, 
L m in = PQi + PQ2 + PQ3, as shown in Fig. [I] 

As a remarkable fact, AE 3 q seems to be relatively well 
expressed as a single-valued function of £ m in- Indeed, 
although there is some visible deviation, the lattice QCD 
data for AE 3 q nearly collapse to a single curve in Fig. [3] 
This is rather nontrivial because A.E3Q depends not only 
on L m in but also on three independent variables. 

For further investigation, we consider the "Mercedes 
form" for the 3Q system as shown in Fig. We define 
Xi = PQ { as the distance between the Fermat point P 
and each quark location Qj. 

After some trials, we finally find that the "inverse Mer- 
cedes Ansatz" defined by the following functional form 
well reproduces the lattice QCD results for the gluonic 
excitation energy A.E3Q in the 3Q system: 

AE 3Q = + (17) 
3 , 

Ly EE ^^i ~ + ^ 0* = P Q.) ( 18 ) 

1=1 

with three parameters, K, G and £. 

Here, we refer to L Y as the "modified Y-type flux-tube 
length", or the "modified Y-length" simply. For £ = 0, 
L Y coincides with the Y-type flux-tube length L m ; n = 

Si=i PQi — Si=i x i- Note that 2Ly expresses the total 
perimeter of the "Mercedes form" shown in Fig. i.e., 

£y = ^E^> ( 19 ) 

since P2Q1 = x\ + £ 2 - 2£xi cos | = x\ — ^x\ + £ 2 etc. 

In Fig.Q we plot the lattice QCD results of the gluonic 
excitation energy AE 3 q ee V 3 q S ' — V3Q ■ against the mod- 
ified Y-length Ly defined in Eq.JTBJ). As a remarkable 
fact, AE 3 q can be plotted as a single- valued function of 
L Y . In Fig. EI w e also plot AE 3 q ee V 3 q S ' — V 3 c q - against 
1 /Ly, since the inverse Mercedes Ansatz corresponds to 
a linear arising behavior in this plot. 
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FIG. 7: The lattice QCD results of the gluonic excitation 
energy A.E3Q = V^q' — V 3 q ' plotted against the modified Y- 
length L— defined in Ea. l|18^ . The dashed curve denotes the 
inverse Mercedes Ansatz with the best-fit parameters listed 
in Table El 



FIG. 8: The lattice QCD results of the gluonic excitation en- 
ergy Ai?3Q = V 3 g Q s ' - V£q plotted against 1/Ly. Tne dashed 
line denotes the inverse Mercedes Ansatz with the best-fit 
parameters listed in Table IV1 



We summarize in Table El the fit analysis with the 
inverse Mercedes Ansatz for the lattice QCD data of the 
gluonic excitation energy A-E3Q at each (3 together with 
the best-fit parameter set, (K, G, £). (In Appendix, 
we summarize various fit analyses for A.E3Q with several 
trial fit functions of {a, b, c} or {xi, X2, X3}.) 



TABLE V: The fit analysis with the "inverse Mercedes 
Ansatz" for the lattice QCD data of the gluonic excitation 
energy A.E3Q at ,3=5.8 and 6.0. The best-fit parameter set 
(K, G, £) is listed in the lattice unit. One finds G(/3 = 5.8) ~ 
0.77 GeV, G(f3 = 6.0) ~ 0.85 GeV, £(/3 = 5.8) ~ 0.116 fm 
and f(/3 = 6.0) ~ 0.103 fm in the physical unit. 



K 



G 



X 2 /N L 



5.8 1.4329(340) 0.5500(105) 0.77 41.2/(24-3)=1.96 
6.0 1.3486(277) 0.4252( 68) 1.03 76.8/(73-3)=!. 10 

In Figs. and |H1 we have added by the dashed line 
the best-fit result of the inverse Mercedes Ansatz with 
the parameters listed in Table Ivl The inverse Mercedes 
Ansatz seems to reproduce fairly well the lattice QCD 
data for the 3Q gluonic excitation energy AE 3 q at least 



2 fm. 

Note again that AE^q is generally a function of the 
three variables {x\, X2, £3}. However, the inverse Mer- 
cedes Ansatz (II Til depends only on Ly, which is a single 
symmetric combination of {xi,X2,xs}. We stress that 
such a simple dependence is rather nontrivial. 

Now, we pay attention to the parameters K, G and £. 
As the physical meaning of the parameters £ and G in the 
inverse Mercedes Ansatz, £ plays the role of an "ultravi- 
olet cutoff" parameter when all of small as Xi < £, 
and the parameter G seems to give an infrared value 
of the gluonic excitation energy. We find a relatively 
good scaling behavior for the parameter set (K, G, £) 
in Table El K(J3 = 5.8) ~ 1.43, K(fi = 6.0) ~ 1.35, 
G(/3 = 5.8) ~ 0.77 GeV, G(f3 = 6.0) ~ 0.85 GeV, 
£(/3 = 5.8) ~ 0.116 fm and £(/3 = 6.0) ~ 0.103 fm in 
the physical unit. 

As a caution, one has to be careful for the argument on 
the infrared behavior of A-E3Q. For, the character of the 
gluonic excitation mode may be changed into the stringy 
behavior in the infrared region, as is actually shown in 
the Q-Q gluonic excitation mode for r > 2 fm. For the 
definite conclusion of the infrared behavior of A-E3Q , we 
need to perform the lattice QCD calculation for AE^q 
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with much larger lattice volume. 

VI. DISCUSSION 

A. Physical implication of the lattice QCD results 

We consider the physical meaning of the present lattice 
QCD result, although the precise physical interpretation 
of the inverse Mercedes Ansatz (II 7|) would be rather dim- 
cult and is an open problem. Since the inverse Mercedes 
Ansatz is described with the modified Y-type flux-tube 
length Ly, the gluonic excitation would be regarded as a 
global excitation of the whole Y-type flux-tube system, 
instead of the partial excitation of each flux-tube as PQ 1; 
PQ 2 or PQ 3 . This would exclude the quasi-fixed edge na- 
ture of the Y-type junction, as was also indicated in sec- 
tiorTvT-A. In fact, the inverse Mercedes Ansatz indicates 
that the gluonic excitation in the 3Q system appears as 
a complicated excitation of the whole 3Q system. 

As a remarkable fact on the absolute value of the glu- 
onic excitation energy, the lowest gluonic-excitation en- 
ergy AE$Q is found to be about 1 GeV or more in the typ- 
ical hadronic scale as 0.5 fm< L m j n < 1.5 fm. In fact, the 
gluonic excitation energy AE^q is rather large in com- 
parison with the low-lying excitation energy of the quark 
origin. (Also for the Q-Q system, a large gluonic exci- 
tation energy is reported in recent lattice studies |17|.') 
Therefore, the contribution of gluonic excitations is con- 
sidered to be negligible and the dominant contribution is 
brought by quark dynamics such as the spin-orbit inter- 
action for low-lying hadrons. 

On the other hand, the gluonic excitation would be 
significant and visible in the highly-excited baryons with 
the excitation energy above 1 GeV. For instance, the low- 
est hybrid baryon |4l| , which is described as qqqG in the 
valence picture, is expected to have a large mass of about 
2 GeV. This lattice QCD result may suggest a large "con- 
stituent gluon mass" of about 1 GeV in the constituent 
gluon picture. 

B. Gluonic excitation and success of quark model 

We consider the connection between QCD and the 
quark model in terms of the gluonic excitation [T5I Il6| . 
While QCD is described with quarks and gluons, the 
simple quark model, which contains only quarks as ex- 
plicit degrees of freedom, successfully describes low-lying 
hadrons, in spite of the absence of gluonic excitation 
modes and the non-relativistic treatment. As for the non- 
relativistic treatment, it is conjectured to be justified by a 
large mass generation of quarks due to dynamical chiral- 
symmetry breaking (DCSB). However, the absence of the 
gluonic excitation modes in low-lying hadron spectra has 
been a puzzle in the hadron physics. 

On this point, we find the gluonic-excitation energy 
to be about 1 GeV or more, which is rather large com- 



pared with the excitation energies of the quark origin, 
and therefore the effect of gluonic excitations is negligi- 
ble and quark degrees of freedom plays the dominant role 
in low-lying hadrons with the excitation energy below 1 
GeV. Thus, the large gluonic-excitation energy of about 
1 GeV gives the physical reason for the invisible gluonic 
excitation in low-lying hadron spectra, which would play 
the key role to the success of the quark model without 
gluonic excitation modes 0, 0] . 

Quantum Chromodynamics 



(Strong) Color Confinement 




(Tightly Stretched) 



Color Flux-Tube Formation 
with a Large String Tension 
a ~ 1 GeV/fm 




Large Excitation Energy 
of the Flux-Tube Vibration 




Large Gluonic Excitation 
Energy AE ~ 1 GeV 




Absence of Gluonic Excitation 
Modes in Low-lying Hadrons 



Dynamical Chiral Symmetry 
Breaking 




Large Constituent Quark Mass 
M 9 ~ 300 MeV 




Non-Relativistic Treatment 
on Quark Dynamics 
in Low-lying Hadrons 



Massive Quark Model for Low-lying Hadrons 

FIG. 9: A possible scenario from QCD to the quark model 
in terms of color confinement and DCSB. DCSB leads to a 
large constituent quark mass of about 300 MeV, which en- 
ables the non-relativistic treatment for quark dynamics ap- 
proximately. Color confinement results in the color flux-tube 
formation among quarks with a large string tension of a ~ 
1 GeV/fm. In the flux-tube picture, the gluonic excitation 
is described as the flux-tube vibration, and the flux-tube vi- 
brational energy is expected to be large, reflecting the large 
string tension. The large gluonic-excitation energy of about 
1 GeV leads to the absence of the gluonic mode in low-lying 
hadrons, which would play the key role to the success of the 
quark model without gluonic excitation modes. 

In Fig. by way of the flux-tube picture, we present a 
possible scenario from QCD to the massive quark model 
in terms of color confinement and DCSB [l(j. On one 
hand, DCSB gives rise of a large constituent quark mass 
of about 300 MeV, which enables the non-relativistic 
treatment for quark dynamics approximately. On the 
other hand, color confinement results in the color flux- 
tube formation among quarks with a large string tension 
of a ~ 1 GeV/fm. In the flux-tube picture, the gluonic 
excitation is described as the flux-tube vibration, and the 
flux-tube vibrational energy is expected to be large, re- 
flecting the large string tension and the resulting short 
string. Due to the large gluonic-excitation energy, which 
is actually estimated as about 1 GeV in lattice QCD, the 
gluonic excitation seems invisible in low- lying hadrons. 
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In this way, the quark model becomes successful even 
without explicit gluonic modes. 



APPENDIX A: FIT ANALYSES FOR THE 
GLUONIC EXCITATION ENERGY 



VII. SUMMARY AND CONCLUSION 

For about 100 different patterns of spatially-fixed 
three-quark (3Q) systems, we have studied the excited- 
state 3Q potential and the gluonic excitation using SU(3) 
lattice QCD with 16 3 x 32 at /3=5.8 and 6.0 at the 
quenched level. We have extracted the excited-state 
potential V 3 q' together with the ground-state potential 
V 3 q S ' by diagonalizing the QCD Hamiltonian in the pres- 
ence of three quarks, for 24 patterns of 3Q systems at 
(3 — 5.8 and for 73 patters at (3 = 6.0. 

We have found that the lowest gluonic excitation en- 
ergy AE 3 q = V£q — V^q' takes a large value of about 
1 GeV at the typical hadronic scale as 0.5 fm < L m ; n < 
1.5 fm. Therefore, we have conjectured that the "hybrid 
baryon" qqqG, which corresponds to the gluonic excita- 
tion mode, appears as the highly-excited baryon with the 
excitation energy of above 1 GeV. 

Next, we have investigated the functional form of the 
gluonic excitation energy A.E3Q in terms of the 3Q loca- 
tion. After some trials with various functions, we have 
found that the lattice data of A^q are relatively well 
reproduced by the "inverse Mercedes Ansatz" , AE 3 q = 
K/Ly + G with the "modified Y-type flux-tube length" 
Ly- This nontrivial behavior of AE 3 q seems to indicate 
that the gluonic-excitation mode is realized as a compli- 
cated bulk excitation of the whole Y-type flux-tube sys- 
tem, instead of the partial excitation of each flux-tube. 

Finally, we have considered the physical consequence of 
the large gluonic-excitation energy, and have presented a 
possible scenario to give a physical reason of the success 
of the quark model for low-lying hadrons even without 
explicit gluonic modes. 

These lattice QCD data of the excited-state potential 
would be useful for the QCD-based construction of the 
refined quark model, which can deal with the gluonic 
excitation modes in hadrons. Our results would be also 
helpful for the comprehension of the nature of the "QCD 
string" . 
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In Appendix, we summarize various fit analyses for 
the gluonic excitation energy in the spatially-fixed 3Q 
system, AE 3 q = V 3 q' — V 3 S q' , in terms of the 3Q spa- 
tial configuration. We examine various trial fit functions 
of {a, b, c} or {xi, x%, x 3 } including two or three free pa- 
rameters, af. Form 1-10 expressed by Eqs.(Al)-(A10). 
({a,b,c} and {xi,X2,x 3 } are defined in sectiodVI) 

Form 1 is the "inverse Mercedes Ansatz" , which is dis- 
cussed in sectiorfvl as the best fit form for the gluonic 
excitation energy Ai?3Q. Form 2 is a simplified form of 
Form 1. We examine Form 3-5 and Form 6-8 as the typ- 
ical single- valued function of L m i n and La , respectively. 
As a possibility, the gluonic excitation energy A£"3Q may 
be controlled by the "typical length" of the static 3Q sys- 
tem, and hence we examine also Form 9 and Form 10, 
regarding max(ii + x 2l x 2 + x 3l x 3 + x\) and max(a, 6, c) 
as the typical length of the Y-type flux-tube system. 

For each of Form 1-10, we show x 2 /^df as the fit 
result for the lattice data of AL3Q together with the best- 
fit parameters in Tables IVTl and IVTfl at /3=5.8 and 6.0, 
respectively. 

As a result, Form 1, the "inverse Mercedes Ansatz", 
is found to be the best fit function for AL3Q, which 
reproduces the lattice data of AS3Q with the smallest 
X 2 /N DF ~ 1. 



TABLE VI: The best-fit analysis for the lattice QCD data 
of the gluonic excitation energy AE^q at /3=5.8 with various 
trial fit functions. For each form, the best-fit parameters at 
are listed with the xV-^df- 



Form 


a% 


<Z2 


0,3 


X 2 /N DF 


(Al) 


1.4329(340) 


0.5500(105) 


0.77 


1.96 


(A2) 


2.0996(501) 


0.4936(117) 


0.55 


2.16 


(A3) 


2.6676(3674) 


1.7127(3568) 


0.4734(339) 


3.44 


(A4) 


-0.2017(2759) 


0.3210(293) 


1.4428(928) 


3.27 


(A5) 


0.3818(827) 


1.4193(1612) 


0.1044(1879) 


3.29 


(A6) 


5.1692(7427) 


3.4586(6932) 


0.4515(369) 


2.18 


(A7) 


-0.0067(5676) 


0.3417(334) 


1.8386(1740) 


2.02 


(A8) 


0.3269(830) 


1.8641(1922) 


-0.0411(2578) 


2.02 


(A9) 


0.9341(213) 


0.6304(89) 




11.71 


(A10) 


0.8065(183) 


0.6352(87) 




12.84 
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• Form 1 (The inverse Mercedes Ansatz) 
ai di 



A-E30 = h a 2 = 



L 



Y 



• Form 2 



A-E3Q — 



Oi 



a 2 



+ a 2 (Al) 



(A2) 
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TABLE VII: Thebest- 
at (3=6.0 with various 



fit analysis for the lattice data of Ai?3Q 
trial fit functions. 



Form 


ai 


a-2 


A3 


X /Ndf 


(Al) 


1.3486(277) 


0.4252(68) 


1.03 


1.10 


(A2) 


1.9810(410) 


0.3860(73) 


1.01 


1.52 


(A3) 


2.3757(2175) 


2.0816(3065) 


0.3794(146) 


2.59 


(A4) 


-0.5066(2125) 


0.2871(164) 


1.0987(459) 


2.53 


(A5) 


0.4284(559) 


1.0722(495) 


0.1760(694) 


2.54 


(A6) 


4.6197(4238) 


4.2849(5709) 


0.3644(156) 


1.57 


(A7) 


-0.3829(4184) 


0.3065(180) 


1.3764(780) 


1.54 


(A8) 


0.3659(538) 


1.3968(512) 


0.0862(919) 


1.54 


(A9) 


0.82575(164) 


0.4846(63) 




8.67 


(A10) 


0.7164(143) 


0.4864(63) 




9.31 



• Form 6 



• Form 7 



• Form 8 



AE- 



3Q 



ai 



La + a 2 



+ a 3 



3Q 



«3 



(L A + ai) Q2 



a 2 



AE 3Q = + as 



(A6) 



(A7) 



(A8) 



• Form 3 



• Form 4 



• Form 5 



AE : 



3Q 



ai 



-^min + ^2 



Ai?3Q = 



a 3 



a 3 



(L min + ai) a 2 



a 2 



AE 3Q - ^-^r + «3 

1^ m i n 



(A3) 



(A4) 



(A5) 



• Form 9 

A-E3Q = 



ai 



max(xi + x 2 ,x 2 + x 3 ,x 3 + xi) 



a 2 (A9) 



• Form 10 



AE : 



3Q 



ai 



max(a, 6, c) 



a 2 



(A10) 
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